Manipulation of photon statistics of highly degenerate chaotic radiation 
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Highly degenerate chaotic radiation has a Gaussian density matrix and a large occupation number 
of modes /. If it is passed through a weakly transmitting barrier, its counting statistics is close to 
Poissonian. We show that a second identical barrier, in series with the first, drastically modifies 
the statistics. The variance of the photocount is increased above the mean by a factor / times a 
numerical coefficient. The photocount distribution reaches a limiting form with a Gaussian body and 
highly asymmetric tails. These are general consequences of the combination of weak transmission 
and multiple scattering. 

PACS numbers: 42.50.Ar, 42.25.Bs, 42.50.Lc 



Chaotic radiation is the name given in quantum optics 
to a gas of photons that has a Gaussian density matrix m . 
The radiation emitted by a black body is a familiar exam- 
ple. The statistics of black-body radiation, as measured 
by a photodetector, is very close to the Poisson statistics 
of a gas of classical independent particles. Deviations 
due to photon bunching exist and were studied long ago 
by Einstein J2|. These are, however, small corrections. 
To see effects of Bose statistics one needs a degenerate 
photon gas, with an occupation number / of the modes 
that is ^ 1. Black-body radiation at optical frequencies 
is non-degenerate to a large degree (/ ~ ^-fi^/kT ^ j^-j^ 
even at temperatures reached on the surface of the Sun. 

The degeneracy is no longer restricted by frequency 
and temperature if the photon gas is brought out of ther- 
mal equilibrium. The coherent radiation from a laser 
would be an extreme example of high degeneracy, but 
the counting statistics is still Poissonian because of the 
special properties of a coherent state m . One way to cre- 
ate non-equilibrium chaotic radiation is spectral filtering 
within the quantum- limited line width of a laser H . This 
will typically be single-mode radiation. For multi-mode 
radiation one can pass black-body radiation through a 
linear amplifier. The amplification might be due to stim- 
ulated emission by an inverted atomic population or to 
stimulated Raman scattering Q. Alternatively, one can 
use the spontaneous emission from an amplifying medium 
that is well below the laser threshold g. 

The purpose of this paper is to show that the statistics 
of degenerate chaotic radiation can be manipulated by 
introducing scatterers, to an extent that would be im- 
possible for both non-degenerate chaotic radiation and 
degenerate coherent radiation. We will illustrate the dif- 
ference by examining in some detail a simple geometry 
consisting of one or two weakly transmitting barriers em- 
bedded in a waveguide (see Fig. n[). For the single barrier 
the photocount distribution is close to Poissonian. The 
mean photocount n is only changed by a factor of two 



upon insertion of the second barrier. But the fluctuations 
around the mean are greatly enhanced, as a result of mul- 
tiple scattering in a region of large occupation number. 
We find that the distribution P{n) for the double-barrier 
geometry is not only much broader than a Poisson dis- 
tribution, it also has a markedly different shape. 

We consider a source of chaotic radiation that is not 
in thermal equilibrium. Chaotic radiation is character- 
ized by a Gaussian density matrix p in the coherent state 
representation ||l|. For a single mode it takes the form 

p= dada*{nfi)^^ cxp{— a* iJ.^^a)\a){a\, (1) 



where /i is a positive real number and \a) is a coherent 
state (eigenstate of the photon annihilation operator a) 
with complex eigenvalue a. If one takes into account 
more modes, a becomes a vector a and p a matrix fi 
in the space of modes. (The factor tt/j, then becomes 
the determinant ||7r/.i||.) We take a waveguide geometry 
and assume that the radiation is restricted to a narrow 
frequency interval Suj around uq. In this case the in- 
dices n,m of a„, p„in label the N propagating waveguide 
modes at frequency loq. 

In thermal equilibrium at temperature T, the covari- 
ance matrix fi = fTL equals the unit matrix 1 times the 
scalar factor / = (g'''^/'^^ — 1)~^, being the Bose-Einstein 
distribution function. Multi-mode chaotic radiation out 
of thermal equilibrium has in general a non-scalar fi. We 
assume that /i. is a property of the amplifying medium, 
independent of the scattering properties of the waveguide 
to which it is coupled. Feedback from the waveguide into 
the amplifier is therefore neglected. 

The radiation is fully absorbed at the other end of the 
waveguide by a photodetector. We seek the probability 
distribution P{n) to count n photons in a time t. It is 
convenient to work with the cumulant generating func- 
tion i^(^) = lnE„e^"P(n)]. For long counting times 



tSio ;:^ 1 it is given by the Glauber formula [filJra 
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lnTr(p:exp[(e«-l)aLtaout]:). (2) 



Here flout is the vector of annihilation operators for the 
modes going out of the waveguide and into the photode- 
tector. The colons : : indicate normal ordering (cre- 
ation operators to the left of annihilation operators). The 
transmission matrix t relates aout = tci to the vector a 
of annihilation operators entering the waveguide. Sub- 
stituting Eq. ([^) for p, we find 

F[£^) — — — In / da.da*\\T:(i\\^^ cxp[—a* fi^^a) 
2n J 
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l)a*tha 



tSuj 



In 1 1 a- (e^- l)ntH\ 



(3) 



In thermal equilibrium, when /x = /I, the determi- 
nant can be evaluated in terms of the eigenvalues r„ of 
the matrix product t^t. The resulting expression 



In thermal equilibrium the deviations from a Poisson dis- 
tribution will be very small, because the Bose-Einstein 
function is <C 1 at optical frequencies for any realistic 
temperature. There is no such restriction on the covari- 
ance matrix /x out of equilibrium. This leads to striking 
deviations from Poisson statistics. 

As a measure for deviations from a Poisson distribu- 
tion we consider the deviations from unity of the Fano 
factor. From Eq. (Eh we derive 
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A Fano factor !F > I indicates photon bunching. For ex- 
ample, for black-body radiation JF = 1-1-/ Q. One might 
surmise that photon bunching is negligible if the wave- 
guide is weakly transmitting, so that N~^Tt tH <C 1. 
That is correct if the weak transmission is due to a sin- 
gle barrier. Then each transmission eigenvalue r„ <ti 1, 
hence J-" « 1. However, if a second identical barrier is 
placed in series with the first one a remarkable increase 
in the Fano factor occurs. 
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has a similar form as the generating function of the elec- 
tronic charge counting distribution at zero temperature 
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where V is the applied voltage. 

The similarities become even more striking if one de- 
rives Eqs. (H) and (||) using the Keldysh approach ||l0| , pj , 
which allows to analyse more general situations. In this 
approach, the cumulant generating function 

Fm^±^Y. f'^ Trln[l 

+ iT„({Gi, e-^«-^/2G,e^«-^/2} + 2)] (6) 

is expressed in terms of the anticommutator of the 
Keldysh matrix Green's functions Gi^r [Q. (The sub- 
scripts l,r refer to particles to the left or to the right of 
the scattering region.) The matrix a^ is the third Pauli 
matrix in Keldysh space. The overall sign is — for bosons 
and -|- for fermions. 

If the eigenvalues of t/xt ' are ^ 1 , we may expand 
the logarithm in Eq. (||) to obtain F(^) = n{e^ — 1), with 
mean photocount n = {tSLu/2TT)Ti fit't. The correspond- 
ing photocount distribution is Poissonian, 
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FIG. 1. Dependence of the Fano factor on the occupation 
number of the modes, for transmission through one (dashed 
line) or two barriers (solid line). The inset shows schemati- 
cally the photodetector (shaded) and the waveguide contain- 
ing one or two barriers. 

Let us first demonstrate this effect for a scalar /i = /I, 
when it has a well-known electronic analogue |]13|,|4| . We 
assume that A^ :^ 1 so that we may replace traces in Eq. 
(^) by integrations over the transmission eigenvalue T 
with density p(T), 



J^ = 1 + / 



J^dTp{T)T^ 
jldTp{T)T 



(9) 



For a single barrier p{T) is sharply peaked at a transmit- 
tance F <C 1. Hence, ^ « 1 for a single barrier. For two 
identical barriers in series the density is bimodal |15[ |, 

p(r) = ^T-3/2(l_r)-l/2^ (10) 

27r 

with a peak near T = and at T = 1. From this distri- 
bution we find that 



r = i + y- 
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While the second barrier reduces the mean photocount 
by only a factor of two, independently of the occupation 
number / of the modes, it can greatly increase the Fano 
factor for large / (see Fig. |l|). From the electronic ana- 
logue (pf) we would find JF = 1 for a single barrier and 
jr=l — i = ifora double barrier [Q. We conclude 
that for electrons the effect of the second barrier on the 
mean current and the Fano factor are comparable (both 
being a factor of two) , while for photons the effect on the 
Fano factor can be orders of magnitude greater than on 
the mean current for / ^ 1. 

The result (^l]) is "universal" in that changing the na- 
ture of the multiple scattering will only change the nu- 
merical coefficient i. For example, multiple scattering 
by disorder would give T = 1 + '^f, in analogy with the 
electronic result ilqJlTIl T = I 
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detector, when 
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For a Lorentzian frequency profile, with maximum /maxj 
one has J^ = 1 + ^/max- 

We now generalize this result to a non-scalar fi. An 
extreme case is a covariance matrix of rank one having 
all eigenvalues ^n equal to zero except a single one. This 
would happen if the waveguide is far removed from the 
source, so that its cross-sectional area A is smaller than 
the coherence area A^ ^- Since Tr (fitH)'^ = (Tr (itH)'^ 
if fj, is of rank one, the Fano factor reduces to J^ = 
1 -f- Tr ntH. The trace of fitH is <C 1 for both a sin- 
gle and double barrier geometry, hence a second barrier 
has no large effect on the noise if A < Ac- 
More generally, for a non-scalar /x the Fano factor (H) 
depends not just on the eigenvalues T„ of tH, but also 
on the eigenvectors. We write tH = U'tU, with U the 
unitary matrix of eigenvectors. We assume strong inter- 
mode scattering by disorder inside the waveguide. The 
resulting U will then be uniformly distributed in the uni- 
tary group, independent of r ||l^. For iV ^ 1 we can 
replace the traces in numerator and denominator in Eq. 
(^ by integrations over U, 
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JdUTi hU^tU 



(13) 



The integrations over U can be carried out easily [Q 
with the result 



.F=1 + (/x)(t) + (/x> 
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(14) 



Here (//P) = iV^Tr fiP, (tP) = iV-^Tr r^ denote the 
spectral moments and {{(J-P)) , ((t^)) the corresponding cu- 
mulants. [For example, ((t^)) = (r^) — (r)^.] 



Instead of Eq. (g_lj) we now have for the double barrier 
geometry a Fano factor 



T= l + i(/x)(l + K), K 
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(15) 



We may estimate the magnitude of the correction k 
by noting that, typically, only Nc ~ A/ Ac eigenvalues 
of fi will be significantly different from 0. If we ig- 
nore the spread among these Nc eigenvalues, we have 
(/x2) « {N/Nc){ny, hence k « T{N/Nc-l). This correc- 
tion will be negligibly small for F ^ 1, unless FA^ > Nc- 
We note that the bimodal distribution dlQ), on which 
the universal result (O) is based, requires the large- A 
regime A^ 3> l/F. For a non-scalar fi the requirement for 
a universal Fano factor is more stringent: 1 ^ FA ^ Nc- 

In the final part of this paper we consider 
the full photocount probability distribution P{n) — 
(27r)^^ /g d^exp[_F(z^) — m^]. For large detection time 
this integral can be done in saddle point approximation. 
The result has the form P{n) = exp[ng{n/n)^ . For small 
relative deviations of n from n the function g(n/n) can 
be expanded to second order in n/n. Thus the body of 
the distribution tends to a Gaussian for t ^ cx), in ac- 
cordance with the central limit theorem. The same holds 
for the Poisson distribution (^. However, the tails of 
P{n) for degenerate radiation remain non-Gaussian and 
different from the tails of Ppoisson(«)- 

Let us first investigate this for a scalar /x = /I. 
Replacing the sum over n in Eq. (0) by the integral 
/„ dTp{T), which is allowed in the large- A^ limit, we find, 
using Eq. (fi(\j, the generating function 
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Nr[l-Jl-iei~l)f]. 



(16) 



The corresponding P{n) is the K-distribution that has 
appeared before in a variety of contexts pp],^. It is 
usually considered only for / <C 1, as is appropriate for 
thermal equilibrium. In the regime 1 <C / ^ n. of interest 
here it has the form 



P(n) = Cexp(-^-^), 



(17) 



with a normalization constant C — {nfn) ^1"^ eKp{2n/ f). 
The essential singularity at n = is cut off below 



n/\/J, where the distribution saturates at P(0) = 
exp(— 2n/-\/7)- One can easily check that the body of 
the distribution (|l^) reduces to a Gaussian with variance 
i/n. 
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FIG. 2. Logarithmic plot of the photocount distribution 
for / = 8 and n ^ oo. The solid curve follows from Eq. (IIm) 
(describing the double barrier geometry) and is very close to 
the large-/ limit (y_7|). The dashed curve is a Gaussian with 
variance (1 + |/)n and the dotted curve is the Poisson distri- 
bution (|7|). (Notice the different vertical scale for the dotted 
curve, chosen such that the Gaussian body of the Poisson 
distribution becomes evident.) 

In Fig. g we compare the distribution ( |l7| ) with a 
Gaussian and with a Poisson distribution, which has 



the asymptotic form Ppo 



(27m) 



-1/2 



exp m 



nln(n/n)]. The logarithmic plot emphasizes the tails, 
which are markedly different. 

For a non-scalar /x we find that the functional form of 
the large-n tail depends only on the largest eigenvalue 
Amax 3> 1 of the Hermitian positive definite matrix tij,t\ 
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(18) 



The number Amax plays the role for a non-scalar fj. of the 
filling factor / in the result (O) for a scalar /x. (For 
broad-band detection Amax should also be maximized 
over frequency.) While the large-n tail is exponential 
under very general conditions, the tail for n <C n has no 
universal form. 

In conclusion, we have calculated the effect of multi- 
ple scattering on the photodetection statistics of radia- 
tion that is both chaotic (like thermal radiation from a 



black body) and highly non-degenerate (like coherent ra- 
diation from a laser). Even for weak transmission there 
appear large deviations of the photocount distribution 
from Poisson statistics, that are absent in the radiation 
from a black body or a laser. They take the form of 
an enhancement of Var n above fi by a factor ex / and 
a slowing down of the large-n decay rate of P{n) by a 
factor 1//. Explicit results have been given for a dou- 
ble barrier geometry, but these findings are generic and 
would apply also, for example, to multiple scattering by 
disorder. Because of this generality we believe that ex- 
perimental observation of our predictions would be both 
significant and feasible. 
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